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ABSTRACT 

The recent detection of polarized radiation from Sgr A* requires a non-thermal electron distribution for 
the emitting plasma. The Faraday rotation measure must be small, placing strong limits on the density 
and magnetic field strength. We show that these constraints rule out advection-dominated accretion 
flow models. We construct a simple two-component model which can reproduce both the radio to mm 
spectrum and the polarization. This model predicts that the polarization should rise to nearly 100% at 
shorter wavelengths. The first component, possibly a black-hole powered jet, is compact, low density, 
and self-absorbed near 1 mm with ordered magnetic field, relativistic Alfven speed, and a non-thermal 
electron distribution. The second component is poorly constrained, but may be a convection-dominated 
accretion flow with M ~ lO _9 M0/yr, in which feedback from accretion onto the black hole suppresses 
the accretion rate at large radii. The black hole shadow should be detectable with sub-mm VLBI. 

Subject headings: accretion, accretion disks — black hole physics — polarization — Galaxy: center 



1. INTRODUCTION 

The nearest supermassive black hole candidate lies at 
the center of the Milky Way galaxy, weighing in at 2.6 x 
10 6 A/ Q , as inferred from motions of stars near the galactic 
center (Ghez et al. 1998; Genzel et al. 1997). The low 
luminosity of the point source associated with the center, 
<; 10 37 erg/s, is a conundrum since accretion from stellar 
winds of neighboring stars should create a luminosity of 
~ 10 41 erg/s. One possibility is that most of the energy is 
carried by the accreting matter into the black hole, as in 
the advection-dominated accretion flow solution (ADAF, 
Narayan & Yi 1994; Narayan, Yi, & Mahadevan 1995). 
Such a situation is achieved when most of the dissipated 
energy is channeled into protons which cannot radiate effi- 
ciently Low efficiency also occurs when gas accretes spher- 
ically and carries its energy in as kinetic energy (Melia 
1992). Alternatively, the accretion rate may be overesti- 
mated, and the emission may be due to a tenuous disk or 
jet (Falcke, Mannheim, & Biermann 1993). 

The radio spectrum of Sgr A* can be described by a 
power law, F v cx v 1 ^ from centimeter to millimeter wave- 
lengths. This is intriguingly close to the spectrum of opti- 
cally thin, mono-encrgctic electrons emitting synchrotron 
radiation (Beckert & Duschl 1997). However, this expla- 
nation is not unique: a self-absorbed source which varies 
in size as a function of frequency may produce a simi- 
lar spectral slope (Melia 1992; Narayan et al. 1995). A 
possible technique to distinguish these models is to mea- 
sure the polarization of the emission: Faraday rotation and 
self-absorption can change the polarization magnitude and 
wavelength dependence (Jones & O'Dell 1977). 

Only recently has linear polarization been detected at 
high frequency by Aitken et al. (2000: A00); previ- 
ous searches at lower frequency showed only upper limits 
(Bower et al. 1999a, b). After correcting for contamination 
by dust and free-free emission, the inferred polarization is 
10-20%, implying a synchrotron origin. This correction 
is made somewhat uncertain by the large beam size. Re- 



markably, the polarization shows a change in position an- 
gle of ~ 90° around 1 mm, which A00 suggest might be 
due to synchrotron self-absorption. We first discuss the 
physics of synchrotron polarization (§2); we then apply it 
to various models in the literature (§3); next we discuss 
a model consistent with all of the observations (§4); and 
finally speculate on the physical implications of this model 
(§5). 

2. SYNCHROTRON THEORY BACKGROUND 

In the synchrotron limit (7 3> 1) for an isotropic 
electron velocity distribution, some analytic results have 
been derived, which we now summarize (Ginzburg & Sy- 
rovatskii 1965 & 1969: GS ). For a uniform slab of elec- 
trons with a power-law distribution, dn e /d , y cx 7"^ (with 
7mm < 7 < Imax such that electrons with "/ min and ^ max 
do not contribute to the frequency of interest), we can re- 
late the magnetic field strength and electron density in 
the slab to the fluid-frame brightness temperature and 
the spectral turnover due to self-absorption. For £ = 2 
and a uniform field B± (projected into the sky plane) 
we find B± ~ 2T n V 2 G and r c ~3x lO^TW^L 
where Tn is the brightness temperature in units of 10 11 
K at the self-absorption frequency v t — 10 12 ^i2 Hz and 
tc is the Compton scattering optical depth of the emis- 
sion region. For v < i/ t , the emission is self-absorbed 
so F v cx z^ 5 / 2 , while above this frequency the emission 
is optically-thin and F v cx j/ 1- ?)/ 2 exp(— vjv max ) where 
v ma x = 3Bj_ejl iax /(4irm e c). 

In the optically-thin regime, the polarization plane is 
perpendicular to the magnetic field with polarization n = 
(£ + l)/(£ + 7/3), up to 100% for £ > 1. In the optically- 
thick regime, n = — 3/(6£+13) (forf > 1/3); the radiation 
polarized perpendicular to the magnetic field is absorbed 
more strongly than the opposite polarization, causing the 
radiation polarized along the magnetic field to dominate, 
switching the polarization angle by 90°, which changes the 
sign of n. Numerical calculations show that the optically- 
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thick polarization peaks at |II| = 20% for £ = 1/3, but 
remains large for < £ < 2. 

To compute the polarization near the self-absorption 
frequency requires a knowledge of the polarized opacity 
and emissivity, A* mi, e_i_ ,||- For £ = 2, these can be ap- 
proximated as (GS): fij_u = (v / Vt)~ 3 (1 ± 3/4) and 

(e±,\\/»±,\\) = 25 t /9(^/^) 5/2 (13 ± 9)/(4 ± 1) where r s 
is the size of the emission region, v t is the frequency for 
which the total source has an optical depth of unity (i.e. 
t = fir s = l/2(^tj_ + (J,\\)r s = 1), St is the source function 
near the frequency v t , and the + or — signs go with the ra- 
diation emitted _L or || to the magnetic held, respectively. 
GS then express the polarization and emission for a slab 
with uniform magnetic field strength and direction, con- 
stant density, and size r s : Ij_ = (ej_//x_i_)(l — exp(— /z_i_r s )), 
7 II = ( e ||/M||)(l-cxp(-A*||r s )),andn= (7 ± -/||)/(/ ± +/||), 
where I±,I\\ are the intensities (erg/cm 2 /s/Hz/sr) with 
polarization perpendicular and parallel to the projected 
direction of the magnetic field on the sky. 

For electron distributions which are highly peaked at 
a single energy (such as mono-energetic or relativistic 
Maxwellian) the polarization for v <J v t is zero. 

The Faraday effect rotates the polarization vector of 
photons emerging from different optical depths by differ- 
ent amounts, causing a cancellation in polarization (Agol 
& Blacs 1996). The differential Faraday rotation angle 



flow of the same luminosity and geometrical thickness. For 
Sgr A*, an accretion rate of ~ 10~( 4 ~ 5 )M Q /yr is inferred 
due to capture of gas in the vicinity of the black hole 
(Quataert, Narayan, & Reid 1999; Coker & Melia 1999), 
which is the value assumed in ADAF models. Assuming 
that the gas falls in at near the free-fall speed, one infers 
an electron density n e — 10 10 cm~ 3 77j_ 5 x~ 3 / 2 and a mag- 
netic field strength of B = W 3 Gm 1 J^x- 5 / i (v A /0.1vff), 
where x is the radius of the emission region in units 



of 



g - GM/c-, m_ 5 
lO~°M /yr, and VA/vff is the ratio of the Alfven speed 
to the free-fall speed. These values imply a total Faraday 
rotation angle at the self- absorption frequency v t of A9 ~ 

10 4 m^ / ' 5 2 ^ 12 2 (w J 4/0.1w//). This value is so large that rota- 
tion of the emitted radiation leads to zero net polarization, 
so ADAFs are in direct conflict with the observed polar- 
ization. Only significant modifications of the model, such 
as a reduction in the accretion rate by a factor of 1CP 3 , 
can reduce the Faraday rotation angle <§C ir. An accretion 
rate of 10~ 8 M Q /yr is consistent with the observed lumi- 
nosity if the accretion flow has a higher efficiency ~2%, no 
longer "advection-dominated." In addition, ADAF models 
assume a Maxwellian electron distribution, which cannot 
produce the observed switch in polarization angle 1 . Fi- 
nally, ADAFs predict a higher self-absorption frequency: 
Ozel et al. (2000) find that u t ~ 5 x lO^m^Hz, which 

implies M ~ 4 x 10~ 7 M Q /yr to be consistent with the 
observed v t ~ 5 x 10 11 Hz. The accretion rate might be 
reduced if there is significant gas lost by a wind or jet 
(Begelman & Blandford 1999; Quataert & Narayan 1999) 
or if the Bondi rate is reduced by heating the infalling gas 
with heat carried outwards by a convection-dominated ac- 
cretion flow, or "CDAF" (Igumcnshchev & Abramowicz 
1999, 2000; Stone, Pringle, & Begelman 1999; Quataert & 
Gruzinov 2000; Narayan, Igumenshchev, & Abramowicz 
2000; Igumenshchev, Abramowicz, & Narayan 2000). 

The model of Melia (1992) is rather similar to the ADAF 
model, and thus suffers the same problems: the high ac- 
cretion rate implies high density which is inconsistent with 
the observed polarization. 

Beckert & Duschl (1997) considered several 1-zone, 
quasi-monoenergetic and thermal emission models for the 
synchrotron emission. These electron distributions do not 
produce a swing in polarization angle by 90 degrees since 
the polarization is suppressed when self-absorbed. Their 
model does produce a self-absorption frequency near the 
correct frequency, however. Falcke, Mannheim, & Bier- 
mann (1993) present a disk-plus-jet model which assumes 
a tangled magnetic field topology which would erase any 
polarization. However, an ordered magnetic field would be 
a small change to their model which might bring it into 
line with the polarization observations. 

4. A PHENOMENOLOGICAL MODEL 

Now, we attempt to construct a model consistent with 
all of the data, using uniform emission regions for simplic- 
ity. Typical optically-thin AGN spectra show £ ~ 2 — 3; 
since £ = 2 is consistent with the polarization from A00, 
we fix £ = 2 in our model fits. The model parameters for 

1 Mahadevan (1999) and Ozel, Psaltis, & Narayan (2000) have added a non-thermal electron component to ADAF models which contributes 
to the flux at wavelengths longer than 2 mm, not at the polarized wavelengths. 



is the accretion rate in units of 



'7„ 



within the source scales as A9 — 3.6 x lQ 2S, T p hotBv~ 
(Jones & O'Dell 1977), where r p hot is the Compton op- 
tical depth of the photosphere. When optically thin, 
Tphot ~ t~c is constant, so rotation is largest at the self- 
absorbed wavelength. When self-absorbed, T p hot of the 
photosphere scales as i^/ 2+2 , so the differential Faraday 
rotation angle oc v^/ 2 (for £ > 1/3), again largest at the 
self-absorption wavelength. The differential rotation at v t 
is A9 ~ 2irg(£ > )(9 b /j mm )t/j mini where 7min is the mini- 
mum electron Doppler factor, g(£) is a dimcnsionless fac- 
tor of order unity, and 9b is the brightness temperature in 
units of m e c /ks- 

3. OBSERVATIONAL CONSTRAINTS ON PUBLISHED 
MODELS 

The observations of polarization in Sgr A* provide the 
following constraints on emission models: 

1) The differential Faraday rotation angle near the self- 
absorbed wavelength must be -C ir. 

2) The electron distribution must be non-thermal since 
the polarization due to a thermal electron distribution is 
suppressed when self-absorbed by a factor of exp(— t). If 
the beam correction by A00 is correct, then II ~ 12% at 
self-absorbed wavelengths, requiring £ <J 2. 

3) The self-absorption frequency must lie near the 
change in polarization angle, ~ 1mm. 

4) The component contributing at lower frequencies 
must have zero linear polarization. 

5) The magnetic field must be ordered to prevent can- 
cellation of polarization. 

These constraints rule out several models proposed in 
the literature, as will be discussed in turn. 

The low efficiency of an ADAF implies a higher accre- 
tion rate and thus higher density than for a high efficiency 
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the polarized component are St = 6 Jy, v t = 550 GHz 
(corresponding to A = 0.55 mm), and v max ~ 5000GHz 
(Figure 1). 

To explain the lack of polarization and spectral slope 
flatter than 5/2, we require an additional component which 
is unpolarized and has a cutoff near f mm so that it 
doesn't dilute the polarization at shorter wavelengths. 
Since Sgr A* has a spectral slope of 1/3 at mm wave- 
lengths and appears to have a spectral turnover at f GHz, 
we model the spectrum as a monoenergetic electron distri- 
bution with energy 7 and zero polarization (due to Fara- 
day depolarization or tangled magnetic field) which be- 
comes self-absorbed at low frequency (Bcckert & Duschl 
1997). For the unpolarized component, we find F v = 
l.3(v/v max ) 1/3 exp(-v/v max )Jy with v max ~ 50GHz, and 
v t ~ 1GHz (Figure 1). 

X(mm) 

1000 100 10 1 0.1 0.01 0.001 




KHz) 



Fig. 1: Polarization and spectral energy distribution of 
Sgr A* compared to model. The dashed line shows the po- 
larized component, the dotted line the unpolarized, mono- 
energetic component, and the solid line the sum of the 
two. The dot-dash line shows the maximum CDAF model 
(assumed to be unpolarized; the total polarization is sim- 
ilar if the CDAF replaces the monoenergetic component). 
The diamonds are the data compiled by Narayan et al. 
(1995), while the asterisks are the data from Bower et al. 
(1999a,b) and A00. 

Figure 1 compares the model to the data. To compare 
the polarization, we have plotted the Stokes' parameter 
that lies at 83°. Remarkably, the polarization should rise 
to ~ 100% at even shorter wavelengths. 

4.1. Physical conditions in Synchrotron emitting regions 

Krichbaum et al. (1998) report a source radius of 55^as 
at 1.4 mm from VLBI observations; this corresponds to 



19r g . The source size may be smaller at higher frequen- 
cies, but we expect the radius of the emission region 
to be greater than the size of the event horizon of the 
black hole, which has an apparent size of ~ 5r g <~ 15^as 
projected on the sky (including gravitational bending, 
Bardcen 1973), so we use an intermediate size in fur- 
ther estimates. The flux of the fitted model at the self- 
absorption frequency, v t = 550 GHz, is ~ 9 Jy. This 
implies a brightness temperature in the emission frame 
T b ~ 1.6 x 10 10 (r s /10r 9 )- 2 r- 1 K, where r s is the size 
of the source (we have assumed the area of the source 
is 7rr 2 ) and T is the bulk Doppler boost parameter. For a 
steeply falling electron number distribution, kT\, ~ 4jm e c 2 
(for £ = 2), where jm e c 2 is the energy of the emitting 
electrons, implying 7 ~ 10(r s /10r g )~ 2 T~ 1 for the elec- 
trons at the self-absorption frequency. Using the for- 
mulae from §2, we find: B ± = 350(r s /10r s ) 4 rG, r c = 
10- 5 (r s /10r 9 )- 8 r- 5 , and lmax _ = 50(r s /10r 3 )- 2 r-\ im- 
plying n e ~ 6 x 10 6 (r s /10r g )~ 9 r~ 5 cm~ 3 . The ratio of 
magnetic to rest-mass energy density is B 2 / (8irn e m p c 2 ) ~ 
l(r s /10r g ) 17 r 9 / 2 for an electron-proton plasma, indicating 
a relativistic Alfven speed. The Faraday rotation angle 
at v t = 5.5 x 10 11 Hz is A0 ~ 350(r s /10r 3 )- 4 r- 2 7 - 3 „, 
assuming B^ ~ B±. For r s ~ 10r g , -f m i n can be as 
large as 4, reducing A9 to 5; for r s ~ 5r g , "f m in can be 
as large as 20 reducing A9 to ~ 0.6. Alternatively, if 
the synchrotron emission is due to a pair plasma, Fara- 
day rotation will be reduced by the ratio of the proton 
number density to the pair number density. The rota- 
tion angle is further reduced at the observed wavelengths 
by a factor ~ v/v t . The high energy cutoff for the elec- 
tron distribution may be due to synchrotron cooling since 
tcooi = 8 x 10 8 7 m ^B- 2 - 6(r s /10r s )- 6 r- 3 sec, similar to 
the dynamical time, to ~ 13x -3 / 2 sec. Given the strong 
scaling of quantities with the unknown r s and T, the above 
estimates can only be improved with future observations. 

The unpolarized emission component dominates at ~ 7 
mm, where Lo et al. (1998) measure a source size of 
~ 5 x 10 13 cm. The self-absorption frequency then requires 
7 ~ 400, B~0.1G, andn e — 4xl0 5 cm~ 3 . Though some- 
what ad-hoc, this model reproduces the spectrum well. 
The Faraday rotation parameter is rather small, so depo- 
larization requires field which is tangled on a scale <~ 100 
times smaller than the size of the emission region. 

4.2. Accretion Component 

We have tried modeling the spectrum of the unpolar- 
ized component with a self-similar, self-absorbed accre- 
tion flow. We used the cyclo-synchrotron emission for- 
mulae from Mahadevan, Narayan, & Yi (1996) and we 
performed the radiation transfer in full general relativity 
(Kurpiewski & Jaroszyfiski 1997). We can place an up- 
per limit on the accretion rate of an ADAF component 
(using the model of Ozel et al. 2000) since its unpolar- 
ized flux must not dilute the polarized component: we 
find Madaf <, 3 x 10" 6 M Q /yr. If the ADAF surrounds 
the polarized emission region, then it will depolarize, so 
the Faraday depolarization places a stronger upper limit 
(§2). We can place a similar limit on the CDAF model 
(using the structure from Quataert & Gruzinov, 2000, 
with equipartition B field and p gas = 2n e kBT e ): we find 
Mcdaf <; 1-5 x 10 _9 M Q /yr; this accretion rate can ac- 
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count for the unpolarized component at v J> 10GHz (see 
Figure 1) and is consistent with the Faraday rotation con- 
straint. The CDAF luminosity is 2 x 10 34 erg/s and the 
self-absorption frequency is ~ 30GHz, so the polarized 
component would be visible through it. Finally, we can 
place a limit on a standard thin disk from the infrared up- 
per limits: we find Mthin ^2x 10~ n M Q /yr; this upper 
limit can be increased to a maximum of 3 x 10~ 7 M Q /yr if 
the inner edge of the disk is truncated at r = 6000r ff . 

5. CONCLUSIONS 

The main success of advection-dominated accretion 
models for Sgr A* is in explaining the high-frequency ra- 
dio spectrum and skirting below the upper limits at in- 
frared frequencies. However, the ADAF model is unpolar- 
ized at the same high frequencies, inconsistent with the re- 
cent detection of linear polarization. We have constructed 
a simple toy model for the millimeter polarization which 
predicts a rise towards shorter wavelengths: polarization 
of - 70% might be seen with SCUBA at 350 fan if this 
model is correct. The lack of polarization and spectral 
slope of 1/3 at wavelengths longer than 2 mm indicates 
that a different physical component may be contributing. 
The presence of two physical components can be confirmed 
by looking for a change in variability amplitude and time- 
scale or source size and morphology around 2 mm. 

The high observed polarization implies a highly ordered 
magnetic field lying near the sky plane. This might be due 
to the poloidal field in a jet (Falcke, Mannheim, & Bicr- 
mann 1993), or due to a toroidal field in a disk component 
seen edge-on. The non-thermal electron distribution might 
be produced by shock acceleration, reconncction, or elec- 
tric field acceleration near the event horizon of a spinning 
black hole (Blandford & Znajek 1977). The Blandford- 
Znajek mechanism can generate a maximum luminosity 
of L BZ ~ 10 37 (£/600G) 2 erg/s (Thorne, Price, & Mac- 
Donald 1986), so the entire polarized luminosity of Sgr A* 
might be powered by black hole spin. 

The dynamics of the emission region will be controlled 
by the ratio of the magnetic field energy density to the 
matter energy density, B 2 / '(8ir pc 2 ); however, this ratio 
scales as r^T 9 / 2 , while T and r s are unknown. Doppler 
boosting decreases the brightness temperature, which re- 
duces Faraday rotation but makes the electrons trans- 
relativistic. Future sub-mm VLBI observations should ac- 
curately measure the r s as a function of frequency, and 
proper motions may constrain T. Also uncertain are the 
pair fraction and minimum electron energy j m i n ■ The pair 
number density can be constrained by measuring the cir- 



cular polarization; without pairs, the circular polarization 
may be as high as a few percent at optically-thin wave- 
lengths (Jones & O'Dell 1977), while pure pair emission 
should have no circular polarization. The pair annihila- 
tion line should be looked for at higher spatial resolution; 
however, it will be strongly broadened by relativistic mo- 
tions of the pairs. Once the source size is known, "f m i n 
and the pair fraction will be constrained by the Faraday 
rotation limit. 

An ADAF model must have a low accretion rate, <J 
lO _8 M0/yr, to be consistent with the lack of Faraday ro- 
tation of the polarized emission. Such a low inferred accre- 
tion rate disagrees with estimates of the Bondi accretion 
rate inferred from stellar winds near the region of the black 
hole. If accretion is episodic due to outer-disk instabilities, 
then the current state might be one of low accretion rate in 
the inner disk. Alternatively, the accretion rate might be 
reduced by depositing energy from the accretion flow in the 
surrounding gas (either through outflow or convection), 
thus increasing the sound speed and decreasing the cap- 
ture rate of gas by the black hole. The accretion flow must 
deposit energy MaGM/ta ~ 6 x 10 35 erg/s, where Ma is 
the stellar mass loss rate which crosses the Bondi radius va 
(Quataert et al. 1999). This can be supplied by accretion 
which releases energy ~ 5 x 10 35 (7y/0.01)m_g erg/s, where 
rj is the efficiency with which accretion deposits energy at 
large radius. As remarked above, a convection-dominated 
accretion flow with M ~ 10~ 9 Af Q /yr can explain part of 
the unpolarized component without diluting the polarized 
emission; the associated convection can carry the required 
energy outward to suppress the Bondi accretion rate. 

Since the self-absorption frequency occurs at ~ 500/im, 
it will be possible to image shadow of a black hole from 
the ground using VLBI, providing a direct confirmation 
of the existence of an event horizon (Falcke, Melia, & 
Agol 2000). Future sub-mm polarimetric VLBI observa- 
tions might show rotation of the polarization angle near 
the black hole, a general relativistic effect which becomes 
stronger for a spinning black hole (Connors, Stark, & Pi- 
ran 1980). 

Eliot Quataert & Andrei Gruzinov have shown me work 
which reaches similar conclusions about ADAFs and the 
need for a very low accretion rate onto Sgr A*. 
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